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1 . 0  INTRODUCTION 


One  of  the  problem  areas  in  mining  and  underground 
construction  is  that  little  is  known  about  the  mass  behavior 
of  rock.  It  is  recognized  that  the  rock  mass  behaves  differently 
from  small  specimens  that  can  be  tested  in  the  laboratory,  but 
the  effect  of  scale  is  not  well  understood,  particularly  under 
generalized  stress  conditions.  Since  large  scale  tests  are 
difficult  and  expensive  even  when  they  are  possible,  investi¬ 
gation  into  the  influence  of  specimen  size  will  help  to  esta¬ 
blish  a  relationship  between  the  behavior  of  small  specimens 
in  routine  laboratory  tests  and  the  mass  behavior  of  rock  in 
the  field.  This  program  vras  intended  to  study  the  scale  effect 
in  triaxial  tests  over  a  range  of  specimen  sizes. 

Two  rock  types,  Charcoal  Black  Granite  (Cold  Spring, 

Minnesota)  and  Indiana  limestone  (Bedford,  Indiana)  were 

1* 

studied  during  this  program.  A  previous  program  considered 
the  granite,  as  well  as  another  rock  type,  so  that  additional 
data  on  large  granite  specimens  was  available.  Four  sizes  of 
specimens,  2-in.  dia,  4-in.  dia,  12-in.  dia,  and  32-in.  dia 
(36-in.  dia.  for  the  limestone)  were  tested  in  triaxial  com¬ 
pression.  All  triaxial  specimens  were  strain  gaged  to  allow 
the  determination  of  elastic  moduli.  The  strain  data  were 
fit  to  models  representing  hydrostatic  and  triaxial  behavior 
to  identify  any  change  in  mechanical  properties  with  size.  The 
failure  data  were  correlated  to  quantify  the  reduction  in 
strength  wij'h  increasing  size. 


Superscripts  indicate  references  listed  at  the  end  of  the 
report 
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2.0  PREVIOUS  STUDIES 
2.1  Effect 'of  Size 

A  number  of  investigators  have  studied  the  strength 
of  pillars  and  coal  cubes  as  a  function  of  size.  Their  results 
are  usually  expressed  as  a  power  function  in  one  of  the 
following  or  similar  forms: 


X 

va 

Rb 

CD 

a*3 

(2) 

oC 

ba 

hP 

(3) 

in  which 

<a  =  failure  stress 
c 

V  =  volume 

R  =  the  pillar  width-height  ratio 
a  =  cube  dimension 

b,  h  =  pillar  width  and  height  respectively. 

2 

Salaman  ar.d  Munro  summarized  reported  values  for  the  exponents 
a,  b,  a  and  (3  in  formulas  of  the  form  of  eq.  (1)  and  (3)  by 
making  use  of  the  relationships  between  (  a,  Q)  and  (a,b) 
implied  by  the  dimensional  terms  V,R,b  and  h.  This  summary 
is  given  as : 


a 

6 

a 

b 

Salaman  and 

Munro z 

-0.66 

±0.16 

0.46 

±0.14 

-0.067 

±0.G48 

0.59 

±0.14 

Greenwald 
et  al3 

-0.183 

0.50 

-0.1-11 

0.72 

Steart^ 

Holland  and 

Gaddy 

-1.00 

0.50 

-0.167 

0.83 
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Salaman  and  Munrc  observe  that  the  variation  in  the  exponents 
may  be  a  function  of  specimen  size,  and  that  there  may  exist  a 
critical  size  above  which  the  effect  of  increasing  volume  is 
negligible. 

Investigators  concerned  with  cube  strength  have  typically 
found  values  of P  in  Eq.  (2)  on  the  order  of  “0.5. 

For  rocks,  the  types  of  relationships  developed  are 

5 

slightly  different.  According  to  the  Weibull  theory  , 

<£  )a  -  Gn) 

m  v 

where  a  =  tensile  or  compressive  strength  of  the  rock  from  a 
standard  laboratory  test, 

«  equivalent  strength  of  the  rock  mass, 

v  =  volume  of  the  rock  mass, 
m 

v  =  volume  of  the  test  sample,  and 

a  =  constant  (with  values  near  10  for  rocks). 

£ 

The  relation  established  by  Protodyakonov  was  of  the  type 

£_  _  i  i  s  (7  -  1) 

a  ^  s  +  a 

m 

where  s  =  spacing  between  major  discontinuj ties  in  the  rock 
mass,  e.g.,  joints,  beds. 

a  =  dimension  of  the  test  specimen,  usually  diameter 
for  compression  tests,  and 
7  =  mass  fract  coefficient  (in  compression:  1-2  for 

igneous  rocks,  1-3  for  competent  sedimentary  rocks, 

3-10  for  weak  rocks;  in  tension:  approximately 
double  these  values). 

7  8  9 

Grobbelaar  based  on  the  work  of  Epstein  ,  Bieniawski  ,  and 

others  found  that  the  formulae  relating  the  modal  strength  of 
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the  weakest  elements  and  its  standard  deviation,  based  on  the 
weakest  link  theory,  are: 

aN  =  au  -  ag  (21ogN)0‘5-  1/2  log(logN)  +  log(4  )  (21ogN)'0’5 

and  a  (N)  =  a  x  ir  (121ogN)”^*^ 

D  O 

or  ag(N)  =  ag(No)  (logho/logN)0*5 

Where  N  =  number  of  flaws  in  the  large  cubic  specimen, 

Nq  =  number  of  flaws  in  the  small  cubic  specimen  or 
unit  cube, 

=  modal  strength  of  the  weakest  link  in  a  sample 
containing  N  elements, 

.a  =  average  strength  of  a  unit  cube  of  material 
containing  Nu  elements , 

ag  =  standard  deviation  of  the  modal  strength  of 
samples  containing  NQ  elements, 

and  o  (N)  =  standard  deviation  of  the  modal  strength 
s 

of  the  weakest  element  in  a  sample  containing 
N  elements. 

These  formulae  are  based  on  the  "weakest  link  theory", 
which  can  be  analyzed  mathematically  if  it  is  assumed  that 
the  frequency  of  occurrence  of  events  is  a  continuous  function 
(e.g.  Weibull^  or  normal  distribution).  The  "links"  in  this 
case  are  the  macroflaws  (or  cracks)  in  the  bulk  material;  not 
the  microflaws. 

Glucklich  and  Cohen  have  indicated  that  effects 

other  than  statistical  exist,  since  the  total  stored  elastic 

energy  increases  with  specimen  volume.  The  energy  released 

at  onset  of  fracture  is  related  to  initiate  fracture;  in  other 

words,  this  reduces  the  rock  strength.  This  phenomenon  has 

12 

been  recently  discussed  by  Baecher 
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2 . 2  Effect  of  Confinement 

There  have  been  numerous  studies  investigating  the 

effects  of  various  aspects  of  confinement  on  rocks.  It  is  not 

intended  to  review  all  of  these  completely  in  this  section. 

Most  of  the  pertinent  work  has  been  briefly  discussed  by 

Swanson^.  The  earliest  experimental  work  was  performed  by 

Adams ^  and  von  Karmen^.  However,  significant  headway  was 

not  made  until  the  initiation  of  work  by  Griggs'^  and  his 
17  18 

coworkers  ' *  .  Since  then,  of  course,  a  number  of  researchers 

have  conducted  various  types  of  studies  under  pressure  several 

19 

of  which  were  presented  at  a  symposium  on  rock  deformation  . 

20 

Baidyuk  has  summarized  some  of  the  Russian  and  American  work. 
Research  in  this  area  is  still  very  active  13,21,22^  All  0f 
this  work  has  been  performed  with  small  rock  specimens,  a  few 
inches  in  diameter.  As  a  result  considerable  light  has  been 
shed  on  the  behavior  of  the  rock  matrix  and  the  criteria  of 
failure.  Refinements  to  the  Griffith  hypothesis  have  been 

no  n /  np 

proposed  *  *  and  appear  to  explain  the  rock  fracture 

process  under  confinement  fairly  well.  The  extrapolation  of 
these  theories  to  larger  rock  masses  is  of  doubtful  value 

26 

and  hence  large  scale  field  testing  has  to  be  resorted  to  . 

The  U.  S.  Bureau  of  Mines  has  undertaken  a  rather  comprehensive 

program  to  collect  field  data  with  the  intention  of  correlating 

27  28 

it  into  a  hypothesis  .  The  contributions  of  Hoek  , 

9  29  30  31 

Bieniawski  ,  Wawersik  ,  Cook  ,  and  Houpert"  to  the  mechanism 

of  brittle  failure  in  rock  deserve  to  be  noted  even  though  the 

studies  were  not  conducted  under  a  confined  state  of  stress. 
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3.0  EXPERIMENTAL  APPARATUS 

In  order  to  conduct  triaxial  tests  for  the  range  of 
specimen  sizes  used  on  this  project,  four  triaxial  cells  were 
set  up  as  shown  below: 


Chamber 

D.  (in.) 

Specimen 
Diameter  (in.) 

Maximum 

Chamber  Pressure 
(ksi) 

Maximum 
Axial  Load 
lb.  xl06 

4.0 

1.95 

30 

0.375 

6.5 

3.65 

30 

0.990 

14.7 

12 

20 

3.40 

48.3 

32  &  36 

20  axial 

10  confining 

36.5 

In  a  standard  triaxial  cell  the  axial  load  is  supplied  by  an 
external  loading  machine.  However,  in  order  to  achieve  the 
large  end*  loads  required  for  the  tests  in  this  program,  these 
chambers  were  separated  into  two  regions  by  sliding  pistons. 
The  general  configuration  is  shown  in  Fig.  1.  One  region 
contained  the  specimen,  and  was  pressurized  to  the  desired 
confining  pressure.  Axial  load  was  transmitted  to  the  rock 
by  the  sliding  piston.  The  maximum  axial  stress  in  the  rock 
depends  upon  the  ratio  of  the  rock  and  piston  areas  and  the 
difference  in  the  confining  pressure  and  the  axial  chamber 
pressure.  The  specimen  stresses  are  given  by  the  following 
equations : 


where 


°2  a3  P3 

ao  -  (prr3)  £ 

+  a  a 

Ap  =  piston  area 
Ar  =  specimen  area 
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P^  =  axial  chamber  pressure 
P^  =  confining  chamber  pressure 

A o  =  deviator  stress 
Op  c*2  and  a ^  =  principal  stresses. 

Reference  to  Figure  1  and  the  above  equations  confirms  that 
if  P^  =  Pg,  A<J  =  0  and  the  specimen  is  under  hydrostatic 
stress  (cr^  =  a2  =  a3^‘  For  ^3  =  sPecimen  i-s  unconfined, 

with  a0  =  a  ’=  0  and  a.  =  A3  =  P,  A  /A. 

1  1  p  r 

3.1  Small  Test  Cells 

The  three  smaller  test  cells  are  incorporated  into 
a  testing  system  with  centralized  controls,  instrumentation, 
and  pumping  systems.  A  schematic  of  this  system  is  shown 
in  Figure  2.  The  tests  conducted  in  this  system  consisted  of 
initial  hydrostatic  loading  up  to  the  desired  confining  pres¬ 
sure,  followed  by  triaxial  compression  at  constant  above 
that  pressure.  At  least  one  load-unload-reload  cycle  was 
observed  for  each  test.  Provision  was  made  for  pressure 
cross-connections  between  the  confining  pressure  and  axial 
pressure  chamber  volumes  to  insure  hydrostatic  conditions 
during  the  hydrostatic  test  phases. 

3 . 2  48-Inch  Test  Cell 

This  test  cell  is  shown  in  Figures  3  and  4.  The  basic 
unit  is  a  48”  I.D.  by  86"  working  length  chamber  having 
20,000  psi  design  working  pressure.  The  chamber  walls  are 
built  up  from  rings  12"  long  which  are  held  in  place  by  a 
3/4"  thick  liner  on  the  inside  diameter.  The  entire  axial 
load  is  carried  by  a  flexible  reaction  frame  which  was  built 
up  from  steel  strap.  The  140  ton  weight  of  chamber  is 
subs tantially  less  than  the  weight  that  would  have  been 
required  by  conventional  chamber  design.  This  is  the 
largest  chamber  currently  available  at  IIT  Research  Institute, 
and  is  capable  of  applying  axial  loads  of  36  million  pounds 
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Fig.  32  IN.  GRANITE  SPECIMEN  BEING 

LOADED  INTO  LARGE  TRIAXIAL  CHAMBER 


Fig  3  32  IN.  GRANITE  SPECIMEN  AFTER 

TESTING  AT  CONFINING  PRESSURE 
OF  2000  PSI 
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to  rock  specimens  as  large  as  3  1/2  ft.  in  diameter. 

As  can  be  seen  in  the  schematic,  the  pumping  and  control 
systems  for  this  unit  are  simpler  than  for  the  small  chambers. 

A  separate  pump  was  used  at  each  end  of  the  chamber.  Accumulators 
were  not  used  because  the  chamber  volume  itself  is  large  in  com¬ 
parison  with  other  available  pressure  chambers. 

The  operation  of  this  chamber  is  similar  to  that  of  the 
smaller  chambers,  except  that  the  turn-around  time  between  tests 
is  on  the  order  of  a  week  instead  of  an  hour.  The  tests  in  this 
chamber  differed  slightly  from  those  in  the  smaller  test  cells. 

In  order  to  maintain  seal  integrity,  a  positive  pressure  differen¬ 
tial  of  at  least  200  psi  was  maintained  across  the  sliding  piston 
during  the  "hydrostatic"  portions  of  the  triaxial  tescs,  and  the 
axial  pressure  was  not  allowed  to  drop  below  400  psi  at  the  bottom 
of  the  load-unload-reload  cycle.  Had  these  precautions  not  been 
taken,  there  would  have  been  danger  of  upsetting  the  piston  seal;, 
thus  aborting  the  remainder  of  the  test. 

3.3  Specimen  Preparation 

Specimens  of  both  rock  types  (Indiana  limestone  and  Char¬ 
coal  Black  granite)  were  delivered  to  IITRI  in  the  form  of  cores 
ranging  from  12  in.  to  36  in.  in  diameter.  Extra  rock  was  included 
to  permit  coring  of  the  2-in,  and  4-in.  dia.  cores.  These  smaller 
cores  were  cut  parallel  to  the  axis  of  the  larger  cores  so  as  not 
to  introduce  complications  because  of  anisotropy.  All  cores  were 
cut  with  a  1:2  dia  to  length  ratio  with  the  exception  of  the 
largest  cores,  which  were  all  60  in.  long.  Thus  the  aspect  ratio 
of  the  Charcoal  Black  granite  cores  was  32/60  =  1:1.87  and  that 
of  the  Indiana  Limestone  was  36/160  =  1:1.67. 

End  preparation  of  the  2-in.  and  4-in.  dia  cores  con¬ 
sisted  of  facing  and  grinding  on  a  lathe  until  the  ends  were 
plane  and  parallel  to  0.0001  inches.  The  larger  cores  were 
capped  in  a  specimen  cage  to  permit  handling.  The  capping 
material  used  was  a  steel-filled  epoxy.  Figure  4A  and  4B 
each  show  an  assembled  32-in.  specimen  in  its  cage.  The  cage 
tie  rods  were  designed  with  end  fittings  that  would  accept 
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tensile  load  only.  Since  the  maximum  tensile  load  that  could  be 
applied  to  the  cage  by  these  tie  rods  corresponded  to  30  psi 
compressive  stress  in  the  rock  specimen,  the  effect  of  the  cage 
on  the  rock  was  negligible. 

An  array  of  foil  strain  gages  were  mounted  on  each  speci- 
men  as  shown  in  Figure  5.  The  number  used  ranged  from  three 
rosettes  on  the  2-in.  cores  to  thirty  rosettes  on  the  32- and  36- 
in.  cores.  These  were  two-element  reosettes  with  1/4  in.  gage 
length  placed  with  the  direction  of  rolling  parallel  to  the 
specimen  axis.  The  gage  placement  procedure  included  the  follow¬ 
ing  s  teps : 

•  grind  the  rock  surface 

#  apply  two  thin  coats  of  gage  cement 

0  visual  inspection  for  voids  in  the  cement  base 

0  affix  gage  and  solder  leads 

e  apply  two  coats  of  gage  coat  for  water  proofing 

0  check  gage  for  continuity  and  response  (soft 
eraser)  and  replace  if  necessary. 

The  instrumented  cores  were  waterproofed  with  latex  ce¬ 
ment  over  the  foil  gages,  and  at  least  two  coats  of  latex  paint 
with  a  thickening  agent  over  the  entire  rock  surface.  Water¬ 
proofing  is  very  important  in  a  tri axial  test  since  both  the 
loading  conditions  and  the  character  of  the  rock  can  be  changed 
by  intrusion  of  fluid  into  the  rock  voids.  In  those  tests  where 
the  specimens  could  not  be  loaded  to  failure  on  the  second  load 
cycle,  the  specimen  was  recovered  intact  and  stripped  of  paint 
for  visual  inspection.  In  several  cases,  traces  of  oil  were 
found  under  the  paint,  but  typically  the  specimens  were  completely 
dry.  In  no  case  was  there  enough  oil  to  do  more  than  dampen  a 
very  small  area  on  the  specimen  surface. 

3.4  Instrumentation  and  Data  Reduction 

The  instrumentation  on  this  program  included  foil  strain 
gages,  and  pressure  gages.  All  data  were  recorded  using  an 
automatic  digital  data  acquisition  system.  The  pressures 
were  read  using  Hiese  bourdon  gages  and  inserted 
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manually  onto  the  data  system  output  as  they  were  read. 

Due  to  the  bulk  of  data  involved,  the  data  were  reduced 
and  plotted  using  an  1108  Univac  computer.  Two  programs  were 
employed  in  the  reduction  of  each  triaxial  test.  The  first  pro¬ 
gram  plotted  the  raw  strain  data  to  show  the  behavior  of  the  in¬ 
dividual  foil  gages.  On  occasion,  the  response  of  individual 
strain  gages  differed  radically  from  the  majority  of  the  gages. 

The  cause  of  this  behavior  was  then  evaluated  as  either  a  key¬ 
punching  error  and  corrected,  or  as  actual  failure  in  the  opera¬ 
tion  of  the  gage,  and  the  data  discarded.  Reasons  for  improper 
gage  operation  include  actual  gage  failure,  or  a  flaw  in  the 
specimen  near  the  gage  location.  After  the  editing  operation, 
the  data  were  reduced  by  a  program  that  solved  the  bridge  equa¬ 
tion  for  each  gage  location  and  printed  out  the  individual  strains 
and  the  strain  averaged  over  the  specimen.  These  data  were  evalu¬ 
ated  for  any  consistent  variation  in  strain  distribution,  such  as 
might  be  caused  by  barreling  of  the  specimen.  If  any  variation 
actually  occurred,  its  magnitude  was  too  small  to  detect  under 
the  normal  gage  to  gage  variation.  In  addition  to  printed  out¬ 
put,  the  following  plots  were  produced  for  each  test: 

Shear  strain  vs.  deviator  stress 
Volumetric  strain  vs.  mean  stress 
Circumferential  and  axial  strain  vs.  axial  stress 
Elastic,  shear  and  bulk  moduli  vs.  mean  stress. 

Data  from  each  test  was  then  assembled  by  hand  to  show  trends 
from  test  to  test. 
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4.0  EXPERIMENTAL  PROGRAM 

The  experimental  work  consisted  of  non-destructive 
tests  and  unconfined  compression  tests  on  2-in,  dia.  cores,  and 
the  major  series  of  triaxial  compression  tests. 

4.1  Non-Destructive  Tests 

A  series  of  non-destructive  tests  were  conducted  on  all 
cores.  These  tests  provide  information  on  the  varibility 
between  specimens  as  well  as  an  indication  of  the  uni  formity 
of  individual  specimens.  These  tests  include: 

•  Schmidt  Hammer 

•  Shore  Sclera cope 

•  Dilatational  Wave  Velocity 


In  general  these  measurements  were  made  on  the  cores  at  the 
locations  of  the  foil  strain  gages.  The  Schmidt  hammer  re¬ 
bound  test  was  found  to  be  strongly  affected  by  the  size  of 
the  small  cores  -  The  hammer  is  in  contact  with  the  rock,  for 
a  relatively  long  time  during  impact,  and  the  rebound  height, 
is  influenced  by  the  method  being  used  to  hold  the  small  cores. 
After  experimenting  with  several  clamping  and  bedding  arrange¬ 
ments,  the  most  consistent  technique  found  was  to  measure  the 
Schmidt  hammer  rebound  on  the  ends  of  the  specimen  rather  than 
the  sides.  The  specimen  being  tested  was  placed  on  a  large 
block  of  the  parent  rock  to  reduce  the  amplitude  of  the  re¬ 
flected  wave.  This  technique  was  used  for  all  2-in.  dia.  and 
4-in.  dia.  specimens  for  the  Schmidt  hammer  tests.  Larger 
cores  were  tested  on  the  cylinderical  surface  at  gage  location. 


The  velocity  of  the  dilatational  wave  was  measured  at 
all  gage  locations  on  all  cores.  The  system  used  for  the 
small  cores  is  shown  in  figure  6A.  The  input  signal  is  a 
single  cycle  square  wave.  This  is  amplified  and  used  to  drive 


a  1  MHz  piezoelectric  transmitting  transducer.  The  signal 
transmitted  through  the  rock  is  received  by  a  similar  transducer 
and  displayed  on  a  dual  beam  oscilloscope  through  an  internal 


variable  delay  line. 


signal  is  compared  with  the  input 
(IT  RESEARCH  INSTITUTE  * 


17 


Fig.  6A  ACOUSTIC  VELOCITY  SYSTEM  FOR  SMALL  CORES 


Fig.  6B  ACOUSTIC  VELOCITY  SYSTEM  FOR  LARGE  CORES 


signal,  and  the  variable  delay  line  adjusted  to  achieve  coin¬ 
cidence.  The  elapsed  time  is  then  read  from  the  delay  line. 

The  system  delay  is  easily  measured  by  removing  the  rock  speci¬ 
men  and  placing  the  transducers  face  to  face.  Figure  6B  shows 
the  system  used  for  the  32-in.  and  36-in.  dia  cores.  A 
mechanical  delay  line  was  introduced  to  allow  the  use  of  rapid 
sweep  rates,  and  a  preamp  boosted  the  received  signal.  The 
operating  frequency  ranged  from  300  kHz  for  the  large  cores 
to  1  MHz  for  the  small  specimens.  A  check  was  made  to  deter¬ 
mine  if  the  frequency  response  of  the  rock  influenced  the 
apparent  time  of  arrival  by  making  measurements  on  a  single 
2-in.  dia  specimen  at  frequencies  between  10  kHz  and  2  MHz. 

The  results  are  plotted  in  Figure  7,  and  indicate  no  variation 
in  time  of  arrival  across  the  entire  frequency  range.  The 
amplitude  of  the  received  signal  is,  of  course,  strongly 
dependent  on  input  frequency. 

4.1.1  Results  of  Non-destructive  Tests -Charcoal  Black  Granite 

Of  the  three  rapid  non-destructive  tests  used  to 
characterize  the  properties  of  the  individual  granite  cores, 
the  sonic  velocity  determinations  were  perhaps  most  useful. 

The  Schmidt  hammer  test  was  found  to  be  size  dependent,  the 
results  of  tests  on  the  2 -in.  dia  cores  being  somewhat  question¬ 
able.  Various  limitations  of  the  Shore  scleroscope  test 
apparatus  make  this  test  inappropriate  for  determining  gross 
properties  in  granite,  unless  test  procedures  are  used  that 
cannot  fit  into  the  category  of  ’’rapid  non-destructive”  test. 

For  best  results,  the  scleroscope  should  be  firmly  held 
vertically  over  the  horizontal  specimen  surface,  as  conditions 
of  plumbness  and  perpendicularity  are  critical.  This  test 
measures  hardness  over  a  very  small  area  on  the  surface  of  the 
specimen,  the  rebound  being  affected  by  the  type  of  individual 
grain  that  is  impacted,  and  the  proximity  of  the  grain  boundary. 
Recommended  procedures  for  good  results  would  involve  cutting 
and  polishing  a  small  specimen  that  can  be  tested  in  a 
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mounted  scleroscope  apparatus.  The  striker  should  be  targeted 
on  individual  grains,  and  a  number  of  replicate  tests  conducted 
to  determine  the  hardness  of  the  individual  constitutent 
minerals  in  the  rock.  This  procedure  was  not  possible  in  this 
program  since  the  individual  cores  could  not  be  sectioned 
prior  to  triaxial  testing.  The  scleroscope  was  hand  held 
vertically  by  reference  to  an  attached  level  bubble,  and  the 
test  cores  rolled  over  the  floor  to  make  the  desired  test 
point  horizontal.  Figure  8  is  a  histogram  showing  the  scleroscope 
readings  for  the  granite.  The  readings  taken  on  12  in.  dia 
cores  are  shaded.  This  size  was  the  only  one  having  two 
distinct  peaks  at  scleroscope  readings  of  50  and  95.  The 
2  in.,  4  in.  and  32  in.  dia.  specimens  displayed  only  the  upper 
peak  at  readings  of  approximately  60  to  90.  Note,  also  that 
nearly  all  readings  above  100  occured  on  the  12-in.  dia  speci¬ 
mens  . 

The  Schmidt  hammer  tests  proved  more  interesting, 
giving  the  first  evidence  that  the  12-in.  dia  granite  had 
properties  significantly  different  from  the  other  sizes.  It 
should  be  remembered  that  the  32-in.  dia  cores  and  the  12-in. 
dia  cores  were  quarried  on  different  occasions,  although  from  the 
same  quarry.  The  2-in.  and  4-in.  dia  cores  were  drilled  from 
a  piece  of  32-in.  dia  stock  during  the  conduct  of  this  research 
program,  and  should  display  properties  similar  to  the  32-in, 
cores.  Figure  9  shows  histograms  of  granite  specimens.  The 
distributions  of  results  on  12-in.  and  32-in.  dia  cores  had 
standard  deviations  near  21/2  Schmidt  hammer  units,  but  the 
means  of  the  two  distributions  differ  by  approximately  two 
standard  deviations,  that  is,  5.1  units.  The  2-in.  and  4-in. 
cores  gave  readings  that  were  influenced  by  their  small  size. 
Various  clamping  and  bedding  arrangements  were  tried  in  an 
effort  to  eliminate  rebound  of  the  cores  and  to  reduce  the 
experimental  scatter  within  each  size.  The  most  consistent 
results  were  obtained  by  testing  the  ends  of  the  cores,  which 
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were  placed  on  heavy  cubes  of  granite.  This  increased  the 
path  length  of  the  stress  wave  in  the  core,  and  reduced  the 
amplitude  of  the  reflected  wave  at  the  risk  of  introducing 
variations  due  to  anisotropy.  The  problem  was  not  completely 
solved,  since  both  2-in.  and  4-in.  dia  cores  displayed  stan¬ 
dard  deviations  of  approximately  4  units,  and  depressed  mean 
values,  especially  in  the  case  of  the  2-in.  dia  specimens. 

The  sonic  velocity  determinations  confirmed  both  of 
the  above  tentative  conclusions,  that  the  12-in.  cores  were 
slightly  more  competent,  and  that  the  variations  in  the  2-in., 
4-in.  and  32-in.  dia.  Schmidt  hammer  readings  were  caused  by 
the  small  specimen  size.  The  first  three  histograms  it- 
Figure  10  show  the  sonic  velocities  measured  in  2-in.,  4-in. 
and  32-in.  dia  cores.  The  measured  velocities  spread  rather 
uniformly,  across  the  range  of  4.3  to  5.1  km/sec.  The 
accuracy  of  the  measuring  system  on  replicate  individual  readings 
was  found  to  be  about  2?0,  so  that  the  observed  variation  is 
inherent  in  the  specimens.  The  last  histogram  in  Figure  10 
shows  all  the  granite  sonic  data,  with  the  normal  distribution 
curve  for  the  12-in.  cores  shown  separately.  Again,  the  12-in. 
dia  cores  are  significantly  more  competent  than  the  other 
granite  cores,  with  sonic  velocities  greater  by  about  207>. 

A  detailed  summary  of  the  mean  sonic  velocities  and 
standard  deviations  for  each  granite  core  is  presented  in 
Table  1.  This  shows  that  much  of  the  scatter  in  Figure  10 
results  from  variation  from  core  to  core.  The  4-in.  and  32-in. 
dia  cores  consistently  have  low  standard  deviations  for  individual 
cores,  so  that  the  normal  distribution  curve  shown  in  Figure  10 
with  standard  deviation  =  0.118  provides  a  good  visualization 
of  the  greatest  possible  measuring  system  errors. 

The  anisotropy  of  the  Charcoal  Black  granite  was 
evaluated  by  checking  the  variation  of  sonic  velocity  on  three 
different  diameters  of  the  2-in.,  4-in.,  and  12-in.  dia  cores. 
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TABLE  1 


SONIC  VELOCITY  FOR  INDIVIDUAL  CHARCOAL  BLACK  GRANITE  CORES 

(km/sec) 


2"  dia 

#6 

#7 

#8 

#9 

#10 

all  2" 

mean 

4.83 

4.79 

4.53 

4.40 

4.75 

4.676 

std.  dev 

- 

- 

* 

- 

... 

0.235 

4"  dia 

#11 

#12 

#13 

#18 

#19 

all  4" 

mean 

4.60 

4.77 

4.78 

4.37 

4.36 

A. 609 

std.  dev 

0.17 

0.19 

0.18 

0.18 

12"  dia 

#45 

#46 

ngU 

#48 

#49 

all  12" 

mean 

5.55 

5.39 

5.27 

5.35 

5.409 

std.  dev 

0.064 

0.055 

Qjfgi 

0.064 

0.041 

_ 

0.118 

32"  dia 

#2 

#3 

#4 

all  32" 

mean 

4.69 

4.97 

4,45 

4.701 

std.  dev 

0.094 

0.085 

0.096 

0.233 
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In  all  cases,  the  variation  between  mean  velocities  on  different 
diameters  was  less  than  the  standard  deviation  of  the  individual 
means.  This  indicates  that  any  consistent  anisotropy  present 
was  smaller  than  the  core-to-core  variations.  The  12~in.  dia 
cores  were  also  checked  to  determine  if  a  consistent  variation 
existed  along  the  length  of  the  cores.  In  this  case  also,  the 
core-to-core  variation  was  large  in  comparison  to  any  consistent 
anistropy.  The  32-in.  dia  cores  were  not  checked  because  the 
appropriate  directions  were  not  marked  on  the  individual  cores 
at  the  quarry.  The  data  is  summarized  in  Table  2. 

4, 1 . 2  Results  of  Nondestructive  Tests -Indiana  Limestone 

The  comments  regarding  the  Shore  scleroscope  and 
Schmidt  hammer  tests  on  granite  apply  to  a  lessor  degree  on  the 
limestone.  The  limestone  does  not  have  the  large  grain  size 
of  the  granite,  so  the  scleroscope  should  show  a  single  peak 
value,  and  the  reduced  sonic  velocity  should  result  in  more 
consistent  rebound  hammer  tests  on  the  small  cores. 


The  results  of  all  Shore  scleroscope  readings  on 
limestone  are  plotted  in  Figure  11.  There  was  little  variation 
related  to  core  size  in  these  tests.  The  tests  scattered 
about  a  mean  value  of  13.95  with  a  standard  deviation  of  4.56 
scleroscope  units. 

The  Schmidt  hammer  results  are  plotted  for  the  individual 
specimen  sizes  in  Figure  12.  The  2-in  and  4-in.  dia  cores 
were  tested  parallel  to  the  axis  in  the  same  manner  as  the 
small  granite  cores,  and  display  little  shift  in  the  mean  values. 

The  sonic  velocity  determinations  are  plotted  for  each 
size  and  for  all  sizes  together  in  Figure  13.  Again,  there  is 
no  significant  variation  between  the  different  sized  specimens. 

The  anisotropy  of  the  Indiana  limestone  was  investi¬ 
gated  by  evaluating  the  three  36-in.  dia  cores.  The  45 
individual  acoustic  velocity  tests  on  these  cores  had  a  mean 
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TABLE  2 


ANISOTROPY  OF  CHARCOAL  BLACK  GRANITE 
(km/sec) 


Radi 

al  Variation 

Axial  Variation 

col.  1 

col.  2 

col.  3 

row  1 

row  3 

row  5 

2"  dia 
mean 

std.  dev 

4.70 

0.327 

4.69 

0.203 

4.64 

0.194 

4"  dia 
mean 
std. dev 

4.61 

0.379 

4.59 

0.139 

4.63 

0.155 

12"  dia 
mean 
std.  dev 

5.42 

0.135 

5.42 

0.113 

5.40 

0.110 

5.39 

0.110 

5 .44 
0.145 

5.42 

0.086 
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Shore  Scleroscope  Reading 


Fig. 11  SHORE  SCLEROSCOPE  READINGS  FOR  INDIANA  LIMESTONE 


_ .  _ _ _  _ _  rnggeasasassiSaE^^^^^B 


value  of  4.021  km/sec,  with  95%  confidence  levels  at  +0.017 
km/sec  and  a  standard  deviation  of  0.058  km/sec.  This  value  of 
4.021  +0.017  km/sec  is  taken  as  the  acoustic  velocity  in 
Indiana  limestone.  These  45  readings  may  be  separated  into  3 
groups  of  15  readings ,  each  group  spaced  at  60-deg  intervals 
around  the  axis  of  the  core,  or  as  5  groups  of  9  readings, 
each  group  spaced  at  10  in.  intervals  down  the  length  of  the 
core.  In  this  fashion  the  variation  in  acoustic  velocity  can 
be  considered  as  a  function  of  direction  in  the  horizontal 
plane,  or  as  a  function  of  the  vertical  depth  in  the  bedding 
plane.  This  nonemclature,  and  the  variation  of  acoustic 
velocities  as  a  function  of  direction  and  depth  are  indicated 
in  Figure  14.  Note  that  in  all  cases,  the  95%  confidence 
levels  for  the  individual  subgroups  include  the  mean  value 
of  4.021  km/sec  computed  for  the  entire  45  reading  group. 

Looking  at  the  variation  of  acoustic  velocity,  it 
appears  that  the  maximum  velocity  would  be  measured  along  a 
diameter  between  2-5  and  3-6,  along  a  quarry  bearing  of  N-80°“W, 
and  the  lowest  velocity  would  lie  along  a  bearing  of  N-10°-E. 

The  maximum  variation  is  on  the  order  of  1-1/2  percent  of  the 
mean  value.  Considering  the  variation  of  acoustic  velocity 
with  depth,  there  may  exist  a  trend  of  increasing  velocity 
with  depth.  The  slight  reversal  of  this  trend  at  row  no.  5 
may  be  a  real  effect  caused  by  the  presence  of  the  bedding  plane 
in  the  formation.  In  any  case,  the  variation  is  on  the  order 
of  1%  or  less.  The  acoustic  velocities  for  all  subgroups  lie 
near  or  within  the  95%  confidence  levels  for  the  entire  group 
of  measurements. 

The  results  of  all  nondestructive  tests  on  both  rock 
types  are  summarized  in  Table  3. 

4.2  Classification  Tests 

In  addition  to  the  nondestructive  tests  to  evaluate 
the  specimen-to-specimen  variability,  petrographic  analysis, 
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Velocity  km/sec  Velocity  km/sec 


Specimen  Orientation 
and  Array  of  Measurement 
Locations 


2 


Direction  of  Measurement  (column  to  column) 


Depth  of  Measurement  (row  number) 

Fig.  14  ACOUSTIC  VELOCITIES  IN  36 -IN.  DIA  INDTANA 
LIMESTONE 


TABLE  3 


SUMMARY  OF  NONDESTRUCTIVE  TEST  RESULTS 
SHOWING  MEANS  AND  STANDARD  DEVIATIONS 


Shore 

sclerescope 

Schmidt 

Hammer 

Sonic 

Velocity 

[km/sec) 

mean 

0 

mean 

0 

mean 

0 

2  in.  dia 
granite 

63.8 

- 

34.9 

4.14 

4.68 

0.24 

4  in.  dia 
granite 

83.1 

- 

47.0 

3.77 

4.61 

0.24 

12  in.  dia 
granite 

75.4 

- 

54.6 

2.65 

5.41 

0.11 

32  in.  dia 
granite 

74.7 

- 

49.5 

2.34 

4.70 

0.23 

2  in. .dia 
limes  tone 

12.7 

3.48 

24.2 

.  _  .. 

2.48 

3.99 

0.10 

4  in.  dia 
limestone 

lif .  6 

4.53 

28.0 

3.55 

4.00 

0.06 

12  in.  dia 
limestone 

14.4 

4,40 

28.2 

3.93 

4.00 

0.06 

36  in.  dia 
limestone 

14.4 

4.99 

29.8 

4.10 

4.02 

0.06 

! 
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tensile  and  unconfined  compression  tests  were  conducted  to 
define  the  rock  types. 

4.2.1  Petrographic  Analyses 


Petrographic  analyses  for  the  two  rock  types  were  con¬ 
ducted  under  the  direction  of  Dr.  A.  Howland  of  Northwestern 
University  in  Evanston,  Illinois.  These  reports  are  reproduced 
below. 

Charcoal  Black  Granite 
Macroscopic 

Gray  medium-grained,  equigranular;  grain  size  generally 
1-3  mm;  light  gray  feldspar,  black  hornblende,  and  dark 
brown  to  black  biotite  recognizable. 

Microscopic 

Testure:  porphyritic:  rectangular  crystals  of 
plagioclase  2, 5-3. 5  rnm.  long  in  a 
subhedral  to  anhedral  interlocking 
aggregate  of  plagioclase,  alkali 
feldspar,  biotite,  hornblende,  and 
quartz  with  grain  size  largely  in 
the  0.3  -  1.0  mm.  range. 

Mineralogy 

Plagioclase  ca.  50%, 

Euhedral  to  subhedral  crystals  in  two  distinct 
size  ranges,  giving  a  porphyritic  texture.  Fine 
polsynthetic  twinning  and  strong  oscillatory 
zoning  ranging  from  about  An25  in  the  cores  to 
An^4  in  the  outer  zones ;  thin  untwinned  albitic 
rims.  Many  inclusions  of  biotite,  hornblende, 
and  magnetite.  Slightly  clouded  with  fine 
sericitic  alteration. 

Alkali  feldspar  ca.  15%, 

Subhedral  to  anhedral  grains,  fine  perthitic 
intergrowths,  grid  microcline  twinning  visible 
in  some  grains . 

Biotite  ca.  15% 

Euhedral  to  subh^^al  plates,  brown,  many  dark 
haloes  around  zircon  inclusion. 

Quartz  ca.  10% 

Strong  tendency  to  interstitial  occurrence, 

I  i T  RESEARCH  INSTITUTE 


35 


filling  in  between  rectangular  feldspar  grains. 
Hornblende  ca.  7% 

Pale  green  to  yellow  green  pleochroism,  com¬ 
plexly  twinned  and  intergrown  with  magnetite, 
some  contain  cores  of  relict  pyroxene  (included 
in  hornblende  percentage) . 

Magnetite  ca.  3 %, 

Separate  grains  and  intergrowths  in  hornblende. 

Apatite 

Many  needles  and  stout  prisms  included  in 
other  minerals. 

Zircon 

Inclusions  in  biotite. 

This  rock  would  classify  as  a  granodiorite. 

Indiana  Limestone 

Macroscopic 

Fine-grained  porous  light  gray  limestone;  abundant  oolites 
and  fossil  fragments  visible  with  the  hand  lens. 

Microscopic 

Clastic  limestone  composed  of  oolites  0.1  -  0.15  mm.  in 
diameter,  fossil  fragments  ranging  up  to  1.5  mm.  in  length 
and  about  0.1  mm.  thick  generally  alined  parallel  to  the 
bedding,  angular  to  subangular  polycrystalline  calcite 
fragments,  and  calcite  cement. 

Oolites  show  both  concentric  and  radiating  structures, 
the  former  much  more  abundant.  Some  have  overgrowths  of 
calcite.  They  form  about  35%  of  the  rock. 

Fossil  fragments  with  either  fibrous  or  granular  internal 
structure  are  about  30%>  by  volume.  Grains  of  polycrystal¬ 
line  calcite  and  clear  single  plates  of  calcity  cement 
are  about  20%,  by  volume. 

Pore  space  is  estimated  as  about  15%,  of  the  thin  section, 
but  some  interstitial  material  may  have  been  lost  in  the 
grinding  of  the  section. 
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4.2.2  Tensile  and  Compressive  Tests 

Tensile  and  unconfined  compression  tests  were  conducted 
on  2~in.  dia  specimens  in  accordance  with  ASTM  specifications. 
The  unconfined  compression  test  specimens  were  ins trumentated, 
and  the  data  included  in  the  body  of  the  analysis  as  triaxial 
tests  having  03=0. 


rock  type 


granite 

limestone 


tensile 

strength 

(psi) 

1340 

290 


compressive 

strength 

(psi) 

18700 

3970 


4.3  Description  of  Constitutive  Equations 

In  order  to  conveniently  describe  and  analyze  the 
rather  large  amount  of  data  generated  during  this  program,  the 
data  for  the  individual  tests  were  fit  to  descriptive  models. 

The  variation  in  the  model  parameters  is  then  analyzed  for  the 
range  of  specimen  sizes  and  confining  pressures  employed  in 
the  study.  The  models  that  were  chosen  are  described  in  greater 
detail  in  a  previous  report,  as  well  as  the  underlying  theory. 
Briefly,  the  constitutive  equations  for  an  isotropic  -inear 
elastic  material  are  written  in  terms  of  the  volumetric  response 
and  the  distortional  response  through  the  use  of  the  deviatoric 
stress  and  strain  tensors,  defined  by: 


s'  . .  =  e.  . 

ij 


t  EKk«ij 


deviator  strain  (1) 


o'  .  .  =  a.  . 

ij 

in  which 


3  akk  6ij 


deviator  stress  (2) 


=  stress  tensor 


=  strain  tensor 


Kronecker  delta 
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The  linear  elastic  isotropic  constitutive  equation  given  by: 


0 .  .  =  A6 .  *  £,  i  +  2  ue.  . 
vij  ij  kk  M  ij 


(3) 


in  which  A  and  j ±  are  Lame's  constants,  are  thereby  transformed 
into 


o .  .  —  k6  .  .  £  .  +  2  G  £  .  . 

ij  ij  tck  ij 


(4) 


in  which  G  and  k  are  respectively  the  shear  and  bulk  moduli, 
given  by: 


k  =  kk 
3e 


kk 


(5) 


„  o'. . 

G  =  J1-! 

2e< . 
J-J 


(6) 


Nonlinear  behavior  is  incorporated  by  using  incremental 
constitutive  equations  and  defining  G  and  k  as  functions  of 
one  or  more  of  the  stress  invariants.  Inelastic  behavior 
may  be  incorporated  by  defining  separate  load  and  unloading 
moduli. 

A  nonlinear  hysteretic  model  used  by  Seaman  and 
Whitman^o  study  the  behavior  of  sand  appears  to  be  suitable 
for  representing  the  hydrostatic  behavior  of  rock.  The 
stress-strain  curve  for  this  model  is  shown  in  Figure  15. 

For  virgin  loading 

°m  "  A1  ev  (7) 

and  for  unloading  and  reloading 
=  A0  (e  -  e  )n 

2  '  tr  ti  h  ' 


m 


v 


'vi- 


CS) 


in  which 


a  is  mean  stress 
m 
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3  seas 


— —  -I  Tr’t.  ifil fr; 


msmsm 


35&Sefc^>s 


ev  is  volumetric  strain 
e  „  is  residual  volumetric  strain 

n 

A^,  A2  and  n  are  material  properties 

This  model  represents  both  the  nonlinear  and  the  inelastic 
behavior  of  the  rock  observed  under  hydrostatic  loading. 

Application  of  the  mathematical  model  in  a  computer  code  would 
perhaps  be  most  convenient  in  an  incremental  form  using  tangent 
values  of  bulk  modulus,  kt.  This  modulus  is  a  function  of 
mean  stress  and  for  a  given  stress  represents  the  slope  of  the 
curve  at  that  stress  level. 

For  loading 

k„  =  d0m  =  n  A}/n  <£n-1)'n  (9) 

t  t"  1  m 

ev 


and,  for  unloading  and  reloading 


,  do 
k  =  _ m 

dev 


=  n  A^n 

z  m 


(10) 


Note  that  in  this  model  the  modulus  for  unload-reload  is 
larger  than  the  modulus  for  virgin  loading  at  any  given 
stress  level.  Only  loading  models  were  generated  in  this  work. 

4.4  Triaxial  Compression 

A  triaxial  compression  test  normally  consists  of  two 

phases.  First,  a  hydrostatic  confining  stress,  is  applied 

to  the  specimen  so  that  the  principal  stresses  are  all  equal 

to  a  .  Then,  two  principal  stresses,  a  and  a.,,  are  kept 
c  2  j 

constant  at  o c  while  the  third  principal  stress,  is 
increased.  For  these  conditions  one  of  the  components  of  the 
deviator  stress  is 


“  oL  -  1/3  (ax  +  2a3)  =  2/3  (a^a.) 


(ID 
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If  the  material  is  isotropic  two  principal  strains,  &2 

and  e^,  are  also  equal,  meaning  that  the  corresponding  compo¬ 
nent  of  deviator  strain  is 

ell  =  £1  "  1/3  (el  +  2e3)  =  2/3  (el  ‘  e3>  (12) 

Thus,  for  triaxial  conditions  the  shear  modulus  may  be 
determined  from 

a'u  _  "  °3  (13) 

2e  *  2(e^  ~  £<j) 


At  low  stress  levels  shear  modulus  increases  relatively 

rapidly  with  mean  stress;  however,  at  higher  stresses  it 

remains  relatively  constant.  Torsional  wave  velocities  were 

measured  in  cylinders  of  different  granites  exposed  to  hydro- 

2  2 

static  pressure  between  1  kg/cm  and  4,000  kg/cm  by  Birch 

O  O 

and  Bancroft!  From  the  wave  velocities  the  modulus  of  rigidity, 
which  is  identical  to  the  shear  modulus,  was  determined  and  is 
shown  as  a  function  of  mean  stress  in  Figure  16.  Also  shown 
in  Figure  16  are  values  of  shear  moduli  for  specimens  of  Char¬ 
coal  Black  granite  determined  from  triaxial  tests  performed 
during  this  study  and  a  previous  study.  The  moduli  determined 
by  Birch  and  Bancroft  and  those  for  the  Charcoal  Black  granite 
are  in  good  agreement  and  increase  approximately  with  the  1/10 
power  of  the  mean  stress. 

Values  of  shear  modulus  shown  in  Figure  16  are  tangent 
values,  Gt,  which  refer  to  the  slope  of  the  deviator  stress 
value.  For  nonlinear  behavior  the  tangent  shear  modulus  is 
taken  as  a  function  of  mean  stress  and  may  be  calculated  from 

G  =  1/2  d^al  "  g3^  (14) 

d(£^  ~  £3) 


The  functional  relationship  with  mean  stress  may  be  written 
as 


Gt  = 


c  a 


0.1 


m 


(15) 
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g ; 


I- 


R 


in  which  c  is  a  constant. 

4.5  Poisson’s  Ratio  and  Young's  Modulus 

Two  constants  are  necessary  to  describe  the  stress - 
strain  behavior  of  an  isotropic  elastic  material.  One  set  of 
constants,  k  and  G,  which  separate  volumetric  and  deviatoric 
behavior  have  been  described  previously.  Another  pair  of 
constants  which  are  commonly  used  are  Poisson’s  ratio,  v,  and 
Young's  modulus,  E.  The  constants  E  and  v  are  derived  directly 

from  a  triaxial  test  in  which  the  lateral  stresses,  and  cj^, 
are  equal  to  zero.  For  this  special  test  condition. 


E  =  1,  and 


v  = 


-  e3 


(16) 

(17) 


For  more  general  test  conditions 

9kG _ _ 

3k  +  G  ,  and 


E  = 


„  3k  -  2G 
V  2 ( 3k+G) 


(18) 

(19) 


The  use  of  tangent  values  of  bulk  and  shear  moduli,  k^  and 


Gfc,  in  equations  (18)  and  (19)  will  result  in  tangent  values 
for  Poisson’s  ratio,  vt,  and  Young's  modulus,  Et. 

For  hydrostatic  loading  the  deviatoric  stresses  and 
strains  are  zero  and  G  is  undetermined.  Therefore,  E  and  v 
are  also  indeterminate  from  a  hydros  tat.  For  a  triaxial  test 
with  tfg  equal  to  zero,  E  and  v  may  be  computed  from  (16)  and 
(17).  For  a  triaxial  test  with  unequal  to  zero  but 
constant,  tangent  values  of  E  and  v  may  be  computed  from 

da, 

(32) 


Et  = 


V 


dp'3 


and 


(33) 
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4.6  Charcoal  Black  Granite 
4.6.1  Hydrostatic  Behavior 

The  hydrostatic  stress -strain  data  for  the  Charcoal 

Black  granite  was  fit  to  the  nonlinear  hysteretic  model  described 

32 

by  Seaman  and  Whitman.  The  model  parameters  were  deter¬ 

mined  by  plotting  log  mean  stress  as  a  function  cf  log  volu¬ 
metric  strain.  Recall  that  this  model  is  given  by 

o  =  A  sn 
m  v 

The  parameters  A  and  n  were  computed  for  each  size  specimen, 
grouping  the  five  tests  at  each  size;  together.  The  values  for 
A  and  n  were  used  to  compute  the  bulk  modulus  by  using 
the  transform 


he 


kt  = 


ra 


»  nA 


h  „  (n-l)/n 


a  £ 


m 


v 


where  kt  is  the  tangent  bulk  modulus  and  the  other  quantities  are 
as  given  above.  The  hydrostatic  stress -strain  behavior  for  the 
four  specimen  sizes  is  shown  in  Figs.  17  to  20,  together  with 
the  models  fit  to  each  specimen  size.  The  data  from  the  2-in. 
dia  specimens  fit  the  model  closely  at  stress  above  about  2  ksi. 
This  deviation  is  expected  because  the  model  predicts  zero 
stiffness  at  zero  pressure.  The  data  from  the  4-in.  dia 
specimens  appeared  to  be  slightly  stiffer  than  the  2-in. 
specimens,  but  the  same  model  parameters  provided  a  reasonable 
fit  to  the  data.  The  12-in,  dia  specimens  were  considerably 
stiffer  than  the  smaller  cores,  and  displayed  less  strain- 
hardening.  The  32-in.  dia  specimens  were  similar  in  response 
to  the  2-and  4-in,  dia  cores.  The  model  for  the  two  smaller 
sizes  is  also  indicated  with  the  mcdel  for  the  32-in.  dia 
cores  in  Figure 20  for  comparison.  The  difference  in  the  two 
models  is  net  large,  and  reference  to  Figure  18  will  show  that 
the  4-in.  dia  data  would  fit.  either  model  equally  well.  The 
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o  =  Aen 
m  v 

□ 

4 
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17  MEAN  STRESS  vs.  VOLUMETRIC  STRAIN  (2"  GRANITE) 
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O  2  ksi 
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Fig.  18  MEAN  STRESS  vs.  VOLUMETRIC  STRAIN  (4M  GRANITE) 
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Fig.  19  MEAN  STRESS  vs.  VOLUMETRIC  STRAIN 
(  12"  GRANITE) 
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Fig.  20  MEAN  STRESS  vs.  VOLUMETRIC  STRAIN 
(32"  GRANITE) 


deviation  of  the  12-in.  cores  from  the  other  three  sizes  may 
be  explained  by  the  fact  that  the  2-  and  4-  in.  dia  specimens 
were  cored  from  a  waste  piece  of  32-in.  dia  core.  The  12-in. 
dia  specimens  were  probably  obtained  from  a  different  part 
of  the  quarry,  and  were  found  to  display  slightly  different 
properties  in  the  non-destructive  tests.  The  fact  of  greater 
stiffness  and  less  strain-hardening  may  indicate  a  relative 
lack  of .macroflaws  in  these  specimens. 

The  tangent  bulk  modulus  as  predicted  by  the  models 
and  the  observed  values  are  shown  in  Figures  21  to  24.  The 
model  predictions  fit  the  data  well  except  at  low  stress 
levels.  There  is  somewhat  more  scatter  in  this  form  of  data 
presentation,  as  may  be  expected  from  the  fact  that  the  values  rep¬ 
resent  differentiation  of  experimental  data  over  short  inter¬ 
vals.  Finally,  the  derived  models  are  compared  in  Figures  25 
and  26.  It  appears  that  there  is  little  effect  of  specimen 
size  on  the  hydrostatic  behavior  of  Charcoal  Black  granite. 

The  greater  stiffness  of  the  12-in.  dia  cores  is  probably 
caused  by  slightly  different  properties  in  the  rock.  This  is 
confirmed  by  the  non-destructive  tests,  which  would  not  be  ex¬ 
pected  to  be  influenced  by  specimen  size. 

4.6.2  Triaxial  Behavior 

The  triaxial  behavior  of  the  charcoal  black  granite 

33 

is  represented  by  a  model  used  by  Birch  and  Bancroft 

Gt  =  co'O'1 

in  which 

Gt  =  tangent  shear  modulus 

a  =  mean  stress 
m 

c  =  model  parameter 

The  triaxial  behavior  for  all  sizes  of  the  granite 
was  represented  by  a  single  model,  the  data  being  shown  in 
Figures  27  to  30.  By  following  the  data  for  an  individual 
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Fig.  28  SHEAR  MODULUS  vs.  MEAN  STRESS  -4  in.  GRAN  IT 


test,  it  can  be  seen  that  the  shear  modulus  does  increase  with 
mean  stress  at  intermediate  stress  levels .  As  yielding  begins, 
the  shear  modulus  decreases  rapidly,  deviating  from  the  model. 
Inspection  of  Figure  27,  in  which  the  effect  shows  most  clearly, 
reveals  that  the  point  at  which  yielding  occurs  is  higher  for 
the  tests  at  higher  confining  stress,  exactly  as  would  be 
expected. 

4.6.3  Young's  Modulus  and  Poisson's  Ratio 

The  values  of  shear  and  hulk  moduli  derived  from  their 
respective  models  may  be  employed  to  predict  values  for  other 
elastic  moduli.  This  was  done  for  Poisson’s  ratio  and  the  elastic 
modulus,  using  the  relationships  given  in  Eq.  (18)  and  (19)  on 
page  43. 

v  -  kG  _  3k~2G 

a  3k+G  and  "  2(3k+G) 

in  which  E  =  Young’s  modulus 
v  =  Poisson's  ratio 
k  =  hulk  modulus 
G  =  shear  modulus 


The  predicted  values  for  the  elastic  modulus  are  shown  in 
Figures31  to  34,  together  with  the  observed  values.  The 
agreement  is  good  at  low  and  intermediate  stress  levels.  At 
higher  stresses.,  yielding  begins,  reducing  the  observed 
stiffness . 

Predicted  and  observed  values  for  Poisson's  ratio  are 
displayed  in  Figuret  35  to  38.  Inspection  reveals  that  the 
trends  predicted  by  the  model  are  present,  but  the  values 
predicted  are  lower  in  some  cases  than  the  observed  values. 

The  model  is  fallacious  at  low  stress  levels,  where  the  low 
predicted  bulk  moduli  act  significantly.  As  yielding  occurs, 
Poisson's  ratio  increases  due  to  dilational  effects,  resulting 
in  values  greater  than  0.5.  Yielding  occurs  at  higher  stress 
levels  on  the  tests  at  higher  confining  pressures.  The 
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Fig.  38^  POISSONS  RATIO  FOR  32"  0  GRANITE 


unconfined  tests  shown  in  Fig.  35  begin  yielding  at  a  mean 
stress  of  2  or  3  ksi  (axial  stress  of  6  or  9  ksi) ,  while  the 
test  at  cr-j  =  10  ksi  did  not  yield  until  it  reached  a  mean 
stress  of  23  ksi  (o;  =  49  ksi)  .  The  larger  specimens  behaved 
in  a  more  or  less  similar  manner,  with  little  influence  of 
specimen  size. 

4.7  Indiana  Limestone 

4.7.1  Hydrostatic  Behavior 

The  limestone  typically  displayed  a  strain-softening 
behavior  under  hydrostatic  loadings  The  mean  stress  vs. 
volumetric  strain  data  were  fit  to  the  nonlinear  hysteretic  model 
with  no  significant  variations  attributable  to  specimen  size. 

The  data  from  most  tests  scatter  about  the  curve  given  by  the 
parameters  below. 

o  =  A  £n 
m  v 

A  =  0 . 490xl06 

n  =  0.75 

Note  that  for  n  less  than  one,  the  model  is  concave  toward  the 
strain  axis.  The  data  for  all  tests  are  shown  in  Figure  39- 
The  model  for  bulk  modulus  derived  for  these  parameters  is 
shown  in  Figure  40  ,  together  with  the  observed  data.  The 
bulk  modulus  is  given  by 

k  =  n  A11’’*'  o  (n~l)/n 
t  m 

in  which  A  and  n  have  the  values  given  above.  This  model  pre¬ 
dicts  a  very  high  modulus  at  low  stress  levels,  which  is  not 
the  case  in  the  real  material.  The  actual  data  appear  to  have 
slight  maxima  in  the  area  of  2  ksi,  and  then  decay  along  the 
curve  given  by  the  model . 

4.7.2  Triaxial  Behavior 

The  behavior  of  the  Indiana  limestone  did  not  correspond 

to  the  one- tenth  mean  stress  model.  The  triaxial  data  for  al1 
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Indiana  limestone  tests  were  plotted  in  the  form: 

AO  =  f  (A  e) 

Aa  =  deviator  stress 

Ae  =  shear  strain 

This  is  shown  in  Figure  41,  and  it  can  be  seen  that  the  data 

6  0  3 

are  bounded  by  the  line  Ao  =  3.55  x  10  (Ae)  °  .  At  very  low 
and  very  high  stress  levels  the  measured  strains  are  larger 
than  would  be  predicted  from  the  model.  The  shear  modulus, 

Gt,  for  all  tests  is  shown  in  Figure  42  as  a  function  of  shear 
strain.  The  actual  moduli  are  nearly  constant  at  shear  strains 
less  than  500  microstrains.  As  the  shear  strain  increases,  the 
shear  modulus  drop  uniformly  until  yielding  occurs .  Above 
the  yeild  point  the  modulus  drops  rapidly.  This  behavior  is 
indicated  by  the  solid  symbols  on  Figure  42,  which  are  the 
trajectory  followed  by  test  38,  a  12-in.  dia  specimen  at  6000 
psi  confining  stress.  The  model  generated  by  Figure  41  can  be 
transformed  to  yield  the  tangent  shear  modulus.  This  is  indi¬ 
cated  on  Figure  42  .  It  can  be  seen  that  the  model  represents 
the  shear  modulus  behavior  of  intermediate  strain  levels. 

The  model  Ao  =  B  (a  e)r  can  also  be  transformed  to 
express  shear  modulus  as  a  function  of  deviator  stress.  This 
transform,  together  with  the  measured  data  are  shown  in 
Figure  43  .  The  data  follow  the  model  at  intermediate  stresses, 
but  roll  off  at  low  stress  levels  and  above  yield. 


4.7.3  Young 8 s  Modulus  and  Poisson's  Ratio 


In  order  to  consider  the  elastic  modulus  and  Poisson's 
ratio,  the  models  for  shear  modulus  and  bulk  modulus  are  com¬ 
bined  to  produce  E  and  v  values  for  comparison  with  the  ob¬ 
served  quantities.  In  the  case  of  the  shear  modulus  model, 
the  deviation  from  the  model  is  extreme  at  low  deviator  stresses. 
The  shear  moduli  represented  by  the  curved  line  in  Figure  43 
would  better  account  for  this  roll-off.  The  models  produced 
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Fig.  43  SHEAR  MODULUS  vs.  SHEAR  STRESS  FOR  ALL  LIMESTONE 


for  elastic  moduli  are  shown  in  Figure  44  as  a  function  of 
mean  stress.  The  dashed  lines  indicate  the  change  in  the  models 
that  result  from  the  use  of  the  curved  line  in  Figure  43  . 

The  data  for  each  level  of  confining  stress  are  compared  with 
these  mc^els  in  Figures  44  through  49  .  As  in  the  case  of  the 
shear  model,  which  established  an  upper  limit  for  the  shear 
moduli,  the  observed  elastic  moduli  are  also  bounded  by  the 
model.  The  material  is  reasonably  well  represented  at  inter¬ 
mediate  stress  levels,  but  is  softer  than  indicated  by  the 
model  at  both  low  and  high  stresses. 

In  like  fashion,  Poisson's  ratio  models  are  shown  in 
Figure  50  and  the  data  at  each  confining  stress  in  Figures  51 
through  55  „  These  data  are  less  consistent,  and  about  all 
that  can  be  said  is  that  both  the  models  and  the  data  increase 
with  mean  stress. 

In  discussing  the  elastic  modulus  and  Poisson's  ratio, 
it  should  be  remembered  that  the  models  for  these  properties 
were  derived  using  the  bulk  modulus  model,  which  was  fit  to 
the  hydrostatic  data  only.  The  models  for  G,  E,  and  v  can  be 
applied  uniquely  only  to  the  triaxial  data. 

4.8  Strength  Properties 

The  strength  of  a  rock  is  typically  displayed  by 
plotting  the  locus  of  failure  points  as  a  function  of  the 
loading  conditions.  The  failure  data  for  Charcoal  Black  granite 
and  Indiana  limestone  are  shown  on  Figures  56  A  and  56  B  as 

f 

functions  of  mean  stress  and  J^. 

4.8.1  Charcoal  Black  Granite 

Sufficient  data  now  exists  to  produce  a  failure  curve 
for  the  small  granite  specimens.  The  2-in.  and  4-in.  dia  cores 
scatter  about  the  curve  indicated  in  Figure  56A  ,  with  no  indi¬ 
cation  of  any  influence  due  to  size  in  this  range.  These  cores 
all  failed  catastrophically,  with  no  indication  of  slippage 
along  joints  or  weakness  planes  prior  to  failure. 
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Fig.  56B  FAILURE  STRESS  FOR  INDIANA  LIMESTONE 


No  failure  curve  is  defined  for  the  32-in.  dia  granite  I 

as  yet.  These  cores  fail  by  a  series  of  slippages,  displaying  1 

ability  to  carry  increased  load  long  after  failure  begins.  In  J 

several  cases,  the  32-in.  dia  cores  survived  two  complete  load-  | 

unload  cycles,  and  were  found  to  be  fractured  only  upon  removal  f 

from  the  test  cell.  The  highest  stress  state  reached  by  these  J 

cores  is  shown  on  Figure  56A  ,  but  failure  is  initiated  at  some  | 

point  lower  on  the  stress  trajectories.  An  evaluation  of  the  j 

scale  effect  can  be  made  by  comparing  the  maximum  axial  stress 
attained  in  the  32-in.  dia  cores  with  the  axial  stress  obtained  I 

by  the  intersection  of  the  stress  trajectories  with  the  failure 
curve  for  small  specimens.  The  following  table  shows  these  results.  < 


Table  4  SCALE  EFFECT  ON  FAILURE  OF  GRANITE 


Failure  data  for 

Failure  curve 

Ratio 

32 -in. dia 

small 

cores 

32-in. 

-r  small  cores 

(.stress  in 

ksi} 

(stress 

in  ksi) 

a 

o 

0- 

0 

a. 

a 

°i 

3 

m 

1 

m 

1 

m 

1 

2 

11.9 

31.6 

15 

41 

0.79 

0.77 

4 

12.5 

29.4 

21.5 

56.5 

0.58 

0.52 

6 

15.0 

33.0 

26 

60 

0.57 

0.55 

It  can  be  seen  that  the  32-in.  dia  specimens  failed  at  one-half 
to  three-quarters  of  the  strength  attained  by  the  2 -in.  and 
4-in.  dia  cores.  The  stress  trajectory  for  the  2  ksi  confining 
pressure  test  is  nearly  parallel  to  the  failure  curve,  and 
one  of  the  weaker  2-in.  dia  specimens  failed  at  nearly  the 
same  point  as  the  32-in.  dia  core.  Tests  were  conducted  on 
32  in.  dia  Charcoal  Black  granite  specimens  on  a  previous 
program  yielding  somewhat  lower  strengths.  These  data  are 
indicated  by  open  symbols  of  Figure  56A,  and  show  the  amount 
of  scatter  observed  in  the  large  specimens.  None  of  the  12-in. 
dia  granite  cores  failed  at  stress  levels  possible  in  that 
test  cell. 
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4.8-2  Indiana  Limestone 


The  strength  of  the  small  Indiana  limestone  cores  ran 
counter  to  the  expected  trend  of  reduced  strength  with  increased 
size.  Inspection  of  Figure 56B  reveals  an  apparent  loss  of 
cohesive  strength  in  the  2  in.  and  4  in.  dia  cores,  as  compared 
with  both  the  12  in.  dia  triaxial  tests  and  the  2  in.  dia 
unconfined  compression  test.  This  may  be  caused  by  distress  (disturbance) 
in  the  limestone  matrix  during  the  hydrostatic  loading  phase 
of  the  triaxial  tests.  Such  an  effect  would  be  greater  for 
tests  run  at  higher  confining  stresses,  and  smaller  for  tests 
run  on  larger  cores,  the  latter  due  to  arching  action  in  the 
distressed  limestone.  That  is,  particles  in  the  outer  layers 
of  the  core  move  into  point-to-point  contact,  so  that  the 
distressed  region  displays  no  change  in  frictional  strength, 
only  slight  volume  change  and  greatly  reduced  cohesion.  In 
this  configuration,  a  portion  of  the  hydrostatic  stress  may 
be  carried  around  the  interior  of  a  large  core  by  frictional 
arching  in  the  outer  distressed  region,  resulting  in  somewhat 
reduced  stresses  and  lesc  cohesion  loss  in  the  interior. 

Whether  this  hypothesis  is  correct  is  not  known  at  this  time. 

An  alternative  cause  for  this  behavior  would  be  the  intrusion 
of  liquid  and  pressure  into  the  pore  volume.  This  is  not  felt- 
to  be  the  cause,  because  the  trend  became  noticeable  early 
in  the  test  series,  and  extreme  precautions  against  liquid 
intrusion  were  taken.  Ail  cores  were  given  two  coats  of  latex 
paint  by  dipping  before  any  triaxial  tests  were  conducted. 

After  the  first  two  triaxial  tests  on  the  small  rocks  were 
conducted,  the  remaining  eight  2  in.  and  4  in.  dia  cores  were 
given  additional  coats  of  latex,  so  that  a  heavy  layer  of 
waterproof  paint  was  built  up  on  the  cores.  Secondly,  had 
pressure  intrudeu  into  the  rock,  the  state  of  stress  wculd  „_ave 
baen  radically  changed,  together  with  the  apparent  mechanical 
response  o f  the  rock.  However,  all  sizes  of  limestone  displayed 
the  same  moduli  and  were  modeled  by  the  same  laws.  Hence  it 
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seems  that  the  only  influence  was  an  actual  loss  of  cohesive 
strength,  probably  due  to  distress  in  the  limestone  matrix 
near  the  surface  of  the  specimens . 

The  36  in.  dia  limestone  cores  failed  at  stress  levels 
ranging  from  85  to  21  percent  of  the  12  in.  dia  core  strengths, 
the  greater  strength  reduction  occuring  at  higher  stress  levels, 
as  was  *:lso  the  case  for  the  charcoal  black  granite. 
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The  limestone  behavior  was  found  to  be  bounded  by  the 
following  model: 

Gt  -  r  B1/W(,"1>/r 

in  which 

=  shear  modulus 

A a  -  deviator  stress 

B  and  r  model  parameters 

The  parameters  for  all  models  are  listed  in  Table  5. 

The  models  were  used  to  compute  values  for  elastic 
modulus  and  Poisson’s  ratio.  The  agreement  with  the  observed 
data  was  good  for  elastic  modulus  at  stress  states  below 
yielding,  and  only  fair  for  Poisson’s  ratio. 

The  only  major  influence  of  specimen  size  on  the  stress- 
strain  behavior  of  either  rock  type  occured  in  the  12-in.  dia 
granite  cores.  These  cores  also  displayed  significantly 
higher  sonic  velocity  and  Schmidt  hammer  readings,  and  the 
difference  in  model  parameters  is  considered  to  be  the  result 
of  slightly  different  rock  properties  rather  than  the  influence 
of  scale. 

5.3  Effect  of  Specimen  Size  on  Strength 

Sufficient  data  is  not  yet  available  to  precisely 
describe  the  scale  effect  on  strength,  or  to  evaluate  the 
statistical  hypotheses  describing  scale  effect.  The  2-in.  and 
4-in.  dia  Charcoal  Black  granite  cores  appear  to  share  a 
common  failure  criteria,  while  the  32 -in.  dia  granite  specimens 
failed  at  stress  levels  only  1/2  to  °  '4  of  the  strengths 
achieved  by  the  small  cores. 

The  36-in.  dia  Indiana  limestone  cores  failed  at 
stress  levels  of  approximately  3/4  those  attained  by  the  12- 
in.  limestone. 
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TABLE  5 

HYDROSTATIC  MODEL  PARAMETERS 


2  in .  and 
4  in.  dia 
granite 

12  in.  dia 
granite 


32  in.  dia 
granite 


all  sizes 
limestone 


A 

n 

2.77 

x  108 

1.73 

4.66 

x  107 

1.30 

4.00 

x  108 

1.73 

0.490 

x  106 

0.75 

model 


a  =  Aen 
m  v 


k„  =  nA”’1  o  <n‘1>/n 
t  m 


TRIAXIAL  MODEL  PARAMETERS 


2  in.  and  4  in. 
dia  granite 

<c  =  1.89xl06 

12  in.  dia 
granite 

c  =  1.895x10 8 

32  in.  dia 
granite 

c  =  1.735xl06 

all  sizes 
limes  tone 

r  =  0.30 

B= 3. 55x10^ 

Model 


Gt  =  c  a 


for  granite 


G^rf1  (Ao)(r-1)/rl 
for  limestone 


5.4  Effect  of  Scale  on  Type  of  Failure 


There  was  a  drastic  change  in  failure  with  increasing 
specimen  size  for  both  rock  types.  The  32-in.  dia  granite  and 
36-in.  dia  limestone  specimens  failed  very  gradually  and 
displayed  increasing  load  carrying  ability  long  after  slippage 
was  initiated.  The  failure  of  the  small  cores  was  very  rapid, 
with  all  strength  lost  immediately  upon  initiation  of  failure. 

If  initiation  of  failure  is  considered  rather  than  maximum 
stress  attained,  the  influence  of  specimen  size  on  strength 
would  be  even  greater  than  that  presented  above. 

5.5  Recommendations 

Additional  data  on  these  and  similar  rock  types  should 
be  accumulated.  The  gradual  failure  and  extended  load  carrying 
ability  of  the  larger  specimens  appears  to  be  similar  to 
failures  observed  in  rock  masses.  It  seems  likely  that  models  for 
large  specimens  may  be  developed  that  will  be  more  effective 
in  predicting  mass  behavior. 

The  small  scale  effect  with  respect  to  elastic  properties 
suggests  that  effort  should  be  diverted  from  strain  gage 
instrumentation  to  obtaining  additional  failure  data.  The 
strain  gages  should  probably  be  replaced  with  deflection  gages 
operating  over  the  length  and  diameter  of  the  specimens,  at 
least  in  the  larger  sizes.  This  will  give  an  indication  of  the 
gross  rock  behavior,  including  slippage. 

The  sonic  velocity  mapping  was  found  to  be  effective 
in  predicting  variations  in  basic  properties,  the  system 
used  being  both  rapid  and  accurate.  The  use  of  this  test 
should  be  increased  in  future  programs .  The  amount  of  Shore 
scleroscope  testing  should  probably  be  reduced  sharply,  and 
not  used  as  an  indicator  of  core  to  core  variability. 
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i 


TEST  NUMBER-  3* 

QIRMETER—  JL2.  M 

CONFINING  PRESSURE _ gflflO  FS1 

INDlflNfVUME_5TUNE_ 


i“ - r - 1 - 1 - 1 - 

.00  -omo  10.00  go  Jio  so.oo  jiojib 

VOUJMETRI C-STRftI N  MICRO-STRflJNSCXlQjM^ 


A-72 


-5  a. 00  -Q-.00  5JUUI  -100..  Ofl.  150,  Ofl  200..00 

VOLUMETRIC-  STRm^MICRO-^TRRI NSC XLd**2A 


A-76 
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TEST  NUMBER  34 

DIRHETEP  12  IN 

+ 

CONFININS  PRESSURE 

+ 

INOIRNR  LIME  STONE 

+ 

+ 

8000  PSI 


- , - | - 1 - 1 - 1 - - 1 

"-12.00  0.00  12.00  24.00  36.00  48.00  60.00 

VOLUMETRIC  STRAIN  MICRO-STRAINS<X10^*3> 


A-78 


DEVIATOR  STRESS  KSI 
i.  oo  o.  ao  o.  oo  r.  oo 


1 


+ 


+ 


« 


+ 


+ 


+ 


TEST  NUMBER  35 
DIAMETER  12  IN 
CONFINING  PRESSURE  10000  PSI 
INDIANA  LIME  STONE 


— - - 1 - 1 - 1 - 1 - 1 

0.00  15.00  30.00  45.00  00.00  75.00 

SHEAR  STRAIN  MICRO-STRAINS<X10«*2> 


:  * 


4-  + 


TEST  NUMBER  35 
DIAMETER  12  IN 
CONFINING  PRESSURE  10000  PSI 
INDIANA  LIKE  STONE 


-15 .00  0.00  15.00  30.00  45.00  60.00  75.00 

VOLUMETRIC  STRAIN  MICRQ-STRflINS<X10**2> 


A- 80 


TEST  NUMBER  39 
DIAMETER-  36  IN 
CONTAINS  PRESSURE  2000  P3X 
INDIANA  LIMESTONE  FAILURE  6335 


>10.00 


0.00  10 .10  29.00  30.00  40.00 

SHEAR  STRAIN  MICRO~STRAINS<X10^2> 


50.00 


A-81 


PEVIflTOR  STRESS  KSI 


DEVinTOR  STRESS  KSI 

.00  .50  1.00  1.50  2.00  2.50  3.00  3.50  1.  00 


